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We extend our previous study of the SU(3) gauge theory with Nf = 2 flavors of fermions in the 
sextet representation of color. Our tool is the Schrodinger functional method. By changing the 
lattice action, we push the bulk transition of the lattice theory to stronger couplings and thus reveal 
the beta function and the mass anomalous dimension j m over a wider range of coupling, out to 
g 2 ~ 11. Our results are consistent with an infrared fixed point, but walking is not ruled out. Our 
main result is that 7 m never exceeds 0.45, making the model unsuitable for walking technicolor. We 
give an extensive discussion of dicretization error. 



PACS numbers: 11.15.Ha, ll.10.Hi, 12.60.Nz 



I. INTRODUCTION 

For some time we have been studying SU(iV) gauge 
theories with fermions in the symmetric two-index rep- 
resentation of color These are among the theories 
that have been proposed [9rfllj as candidate models for 
walking technicolor [12Hl5l | . Here we present our most re- 
cent work on the SU(3) gauge theory with Nt = 2 flavors 
of fermions in the sextet representation. 

A technicolor theory must supply Goldstone bosons in 
order to generate masses for the weak vector bosons. To 
that end, the theory must break chiral symmetry sponta- 
neously. In order to generate quark and lepton masses as 
well while avoiding large effects of flavor-changing neutral 
currents, one demands the additional property of walk- 
ing. Here a large separation between the technicolor scale 
Arc and the "extended" technicolor scale Aetc is at- 
tained by having a near-zero of the beta function; the 
running coupling stalls for many decades in energy be- 
fore chiral symmetry breaking sets in at large distances. 
Furthermore, the effect of the large ratio Aetc / Arc on 
the technifermion condensate has to be enhanced by a 
large anomalous dimension "f m of the mass operator tpt/j. 
Current estimates [16| require 7 m ~ 1. In fact 7,„ = 1 
emerges from analyses based on the gap equation for 
walking technicolor [H, [3| (see also [13] )■ 

Since the beta function and 7 m are the important in- 
gredients of walking technicolor, our work has focused 
on measuring them. We do this using Schrodinger- 
functional (SF) techniques [l8l - |22l | , which are a lattice im- 
plementation of the background-field method. For other 
instances of the Schro ding er functional applied to tech- 
nicolor candidates, see |23l - [29| . 

Our first effort Q used Wilson's fermion action with 
an added clover term [3(| to reduce 0(a) effects, and 
was limited to lattices of linear size L = 4a, 6a, and 8a. 



The result was a discrete beta function that appeared 
to cross zero at a renormalized coupling g 2 ~ 2.0, in- 
dicating an infrared fixed point (IRFP). An IRFP indi- 
cates conformal physics at large distances, the antithe- 
sis of confinement. Intending to understand (and re- 
duce) discretization effects, we then 0] went to larger 
lattices, L/a — 6, 8, 12 lfy and began using hypercu- 
bic smearing — fat links [3lL l32j — in the fermion action. 
This work showed that the IRFP of Ref. [l[ was but a 
lattice artifact. Thanks to its numerical stability, the 
fat-link action also enabled us to simulate at stronger 
couplings, out to g 2 ~ 4.6, corresponding to a bare cou- 
pling ft — 4.4. At stronger bare couplings we encountered 
a phase transition that makes it impossible to tune the 
hopping parameter k so as to make the quark mass zero. 
The result of Ref. [H, then, is a beta function that is 
smaller in magnitude than the two-loop result but that 
does not cross zero in the accessible range of couplings. 

In Ref. [H we also calculated the anomalous dimension 
7 m according to the method of [III l33T - [35l | . We found 
that 7 m first follows the one-loop curve but its rise slows 
at strong couplings so that 7, n < 0.6 in the range of 
couplings that we could reach. 

In our work on the SU(4) gauge theory with two-index 
(decuplet) fermions @,H|, we encountered the same phase 
transition that prevents simulation of the massless the- 
ory. We found that augmenting the pure-gauge part of 
the action with a new fat-plaquette term can move this 
transition to stronger bare coupling and stronger renor- 
malized coupling as well. In the present paper, we present 
the result of applying this strategy to the SU(3) theory. 
The new action enables us to reach g 2 ~ 11, which is 
in the vicinity of the zero of the two-loop beta function, 
discussed by Caswell (36| and by Banks and Zaks |37| . 

While at first glance our calculated beta function 
crosses zero near the two-loop zero, further analysis 
shows that this result is not stable under extrapolation 
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to the continuum limit. It is quite possible that the beta 
function of this theory approaches zero in the range of 
coupling that we can study, and then veers away from 
zero at yet stronger couplings — much as envisioned in 
scenarios of walking. 

Whether the beta function crosses zero or no, we find, 
in agreement with our earlier work, that the anomalous 
dimension has left the one-loop curve and leveled off. The 
limit we set at our strongest coupling is 7,„ < 0.45. Our 
earlier work, based on the pure plaquette gauge action, 
quoted values approaching 0.6, but this was without at- 
tempting to remove lattice artifacts. We report here an 
analysis of lattice artifacts that allows smooth extrapola- 
tion to the continuum limit. The results for both actions 
move downwards and the error bars grow under extrapo- 
lation. As a result, the results for the two actions are not 
in violent disagreement. The small value of j m near the 
(real or approximate) fixed point spells trouble for any 
use of the present theory as walking technicolor. 

The plan of this paper is as follows. In Sec. [TT] we 
present the improved lattice action, which is the only dif- 
ference in the simulation method between this paper and 
Ref. [i| • For other details of our simulations we refer the 
reader to Ref. Q. We also briefly discuss the ensembles 
we generated, and how we deal with the autocorrelations 
of our observables. We proceed in Sec. IIIII to present 
our results for the running coupling and the beta func- 
tion. We re-analyze the data of Ref. [3| according to the 
lights of our later paper on the SU(2) theory @, where 
we learned to take advantage of the slow running in order 
to make maximum use of the several lattice sizes in play. 
Naturally, we add in the results of the new simulations 
obtained with the augmented gauge action. Our first re- 
sults show an apparent zero in the beta function; this 
will not survive our analysis of discretization error later 
in the paper. Section IIVI contains our numerical results 
for the mass anomalous dimension j m . We study the 
finite-lattice corrections to both j m and the beta func- 
tion in Sec. [V] The former reaches a smooth continuum 
limit, albeit with larger error bars than those of the raw 
results of Sec. IIVI the latter is not so well-behaved. In 
Sec. I VII we summarize our results and place them in the 
context of other work. In the appendix we examine the 
fermion contribution to the one-loop SF coupling and ex- 
plain that it does not provide any guidance for the range 
of couplings and volumes that we study. 



II. LATTICE ACTION AND ENSEMBLES 

Our action contains a fermion term and two pure gauge 
terms. The fermion term ipDptp is the conventional Wil- 
son action, supplemented by a clover term [301 ] with co- 
efficient csw = 1 [Hi • The gauge links in the fermion ac- 
tion are fat link variables V^(x). The fat links are the nor- 
malized hypercubic (nHYP) links of Ref. 31], where for 
each link (x,/j,) one takes the weighted average V fl (x) of 
links in neighboring hypercubes with weights ai = 0.75, 



a2 = 0.6, «3 = 0.3, reunitarized and subsequently pro- 
moted to the sextet representation. 
The gauge action is 
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J2 Re Tr V„ (x) V v (x + ft) Vj[ (x + v)V} (a 



(1) 



wherein the first term is the usual sum of plaquettes of 
fundamental thin link variables, while the second term 
contains plaquettes made of fat links in the sextet repre- 
sentation as in the fermion action. N — 3 is the number 
of colors while df = 6 is the dimension of the fermion 
representation, here the sextet. The weak-coupling ex- 
pansion of Sq gives the effective bare coupling [7] , 



1 

71 



2N 



T, 



6P6 



6 



(2) 



where Tq — 5/2 is the group trace in the fermions' rep- 
resentation. 1 

As before, we employ the hybrid Monte Carlo (HMC) 
algorithm in our simulations. The molecular dynam- 
ics integration is accelerated with an additional h eavy 
pseudo- fermion field as suggested by Hasenbusch |39| . 
multiple time scales and a second-order Omelyan 
integrator [4l|- We determined the critical hopping pa- 
rameter k c — k c (P) by setting to zero the quark mass, as 
defined by the unimproved axial Ward identity on lattices 
of size L — 12a. 

Without a systematic search, we found that choosing 
/3e = +0.5 removes the strong-coupling phase transition 
so that we can run at bare couplings down to j3 = 2.0; at 
smaller /3 the acceptance deteriorates rapidly, especially 
for larger volumes, so that we did not go far enough to 
find out if and where the strong-coupling transition turns 
up. At f3 = 2.0, the running coupling for L = 6 turns 
out to be g 2 ~ 11. This is close to the two-loop Banks- 
Zaks zero at g 2 = ^(16tt 2 ) ~ 10.6. We list in Table [Q 
the values of ((3,k c ) and the number of trajectories run 
at each volume, along with the length of the trajectories 
and the acceptance. Poor acceptance forced us to shorten 
the trajectory length in many cases from the usual value 
ofl. 

The observables we measure are the (inverse) SF run- 
ning coupling, 1/g 2 , and the pseudoscalar renormaliza- 
tion factor, Zp. (We measure Zp on the same config- 
urations used to determine 1/g 2 .) Both of them turn 
out to have long autocorrelation times. We monitored 
and controlled this problem by running 4 or 8 streams in 
parallel at each j3 and L. After analyzing each stream 



We presented a test of weak-coupling universality in our paper 
on the SU(4) theory p|. 
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TABLE I: pe = 0.5 ensembles generated at the bare couplings 
(P, k c ) for the lattice sizes L used in this study. Listed are 
the total number of trajectories for all streams at given (P, L), 
the trajectory length, and the HMC acceptance. 



a 
P 




T In 


trajectories 

1 t~ n nn qs} Ti H q 1 


trajectory 

1 PTi cH~ n 


acceptance 


3.5 


13349 


6 


74.8 


1.0 


0.93 






8 


15.5 


0.5 


0.97 






10 


8.0 


1.0 


0.75 






12 


37.0 


0.5 


0.88 


2.5 


0.13991 


6 


8.8 


1.0 


0.61 






8 


14.3 


1.0 


0.43 






10 


20.4 


0.5 


0.75 






12 


35.6 


0.5 


0.50 






14 


11.6 


0.5 


0.61 






16 


17.1 


0.5 


0.48 


2.0 


0.14273 


6 


17.2 


1.0 


0.65 






8 


14.2 


0.5 


0.61 






10 


12.0 


0.5 


0.67 






12 


13.4 


0.5 


0.48 






16 


29.6 


0.4 


0.38 



TABLE II: /3(j — ensembles generated at the new lattice 
sizes L = 10a and 14a, which we add to the ensembles listed 
m Ref. Q . Columns as in Table [I] 
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FIG. 1: Running coupling 1/g 2 vs. a/L. The crosses are 
from simulations with /3e = (Ref. [4| and Table HV)l : top to 
bottom, P = 5.8, 5.4, 5.0, 4.8, 4.6, and 4.4. The circles are 
from simulations with f3& = 0.5 (Table [TTlj) : top to bottom, 
P = 3.5, 2.5, and 2.0. The straight lines are linear fits [Eq. ©] 
to each set of points at given (P,Pe); the slope gives the beta 
function. The dotted line shows the expected slope from one- 
loop running. 



p 




L/a 


trajectories 


trajectory 


acceptance 








(thousands) 


length 




4.8 


0.13173 


10 


9 


1.0 


0.95 


4.4 


0.13510 


10 


16 


1.0 


0.80 






14 


12 


0.5 


0.85 



separately, we fit the results of the streams together to a 
constant. We demanded that the % 2 /dof of the constant 
fit not exceed 6/3 for 4 streams, or 10/7 for 8 streams. 
For the largest volume L = 16a at the strongest coupling 
P = 2.0, we were not able to overcome the autocorrela- 
tions in 1/g 2 even with nearly 30,000 trajectories. We 
therefore omit this point from the analysis of the run- 
ning coupling. The autocorrelations in Zp, on the other 
hand, did allow a consistent determination, and thus we 
keep this point in the analysis of the mass anomalous 
dimension. 

In this paper we present our new results, obtained with 
Pa = 0.5, alongside the /3 6 — results presented in our 
earlier paper |4|. For the more extensive study of dis- 
cretization error in Sec. El we supplemented the data of 
Ref. [4| with simulations on new lattice sizes L = 10a, 14a 
at two values of P in strong coupling, as shown in Table [TT1 



III. BETA FUNCTION 

The computation of the running coupling proceeds ex- 
actly as described in Ref. with the same boundary 
conditions on the fermion and gauge fields. In brief, one 
imposes Dirichlet boundary conditions at the time slices 
t = 0, L, and measures the response of the quantum effec- 
tive action. The coupling emerges from a measurement of 
the derivative of the action with respect to a parameter 
T\ in the boundary gauge field, 

/95g ( 1 d(D F D F ) 1 \\ 
9 2 (L) \d V T \D F dri D F )/ n=o 

(3) 

The constant K can be calculated directly from the clas- 
sical continuum action. Only go, Eq. @, appears in the 
latter, which assures that K = 12ir regardless of Pe- 

We presented in Ref. [H the values of the running cou- 
pling g 2 for a number of values of (P, k c ) with /3g = 0. 
Our new results for Pq = 0.5 are shown in Table Hill with 
supplementary data on new volumes for p$ — shown in 
Table HVl Both sets are plotted in Fig.Q] 2 



2 We have dropped from consideration the data given in Q] for 
(fS, Pe) = (4.3, 0) since these were taken in a metastable state 
beyond the first-order boundary. 
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TABLE III: Running coupling, Eq. Q, evaluated at the bare coupling (/3, k c ) with /3g = 0.5 on lattices of size L. The omission 
of the result for L = 16a at /3 = 2.0 is explained in the text. 



p 








i/g 2 








L = 6a 


L = 8a 


L = 10a 


L = 12a 


L = 14a 


L = 16a 


3.5 
2.5 
2.0 


0.2918(13) 
0.1454(32) 
0.0915(32) 


0.2859(48) 
0.1433(45) 
0.1023(37) 


0.2888(53) 
0.1424(42) 
0.0942(56) 


0.2703(58) 
0.1517(36) 
0.1057(51) 


0.1647(65) 


0.1449(68) 



TABLE IV: Running coupling for the strong-couplingcases with (5q = 0. Data for L = 10a, 14a are from the new simulations 
listed in Table [Tit while the other data are from Ref. [j. Data for j3 = 4.6 were taken with k — k c = 0.13320. 



p 








1/3 2 






L = 6a 


L = 8a 


L = 10a 


L = 12a 


L = 14a 


L = 16a 


4.8 


0.3151(50) 


0.3002(28) 


0.3013(61) 


0.3008(45) 






4.6 


0.2692(27) 


0.2615(41) 




0.2525(50) 




0.2289(56) 


4.4 


0.2191(24) 


0.2103(34) 


0.2021(66) 


0.2080(50) 


0.2235(105) 


0.2119(53) 



It is convenient to define the beta function /3(u) for 
u = 1 jg 2 as 

0(u) = = 2/3(.g 2 )/.g 4 = 2 U 2 /3(1/ U ) (4) 

aiogL 

in terms of the conventional beta function (3(g 2 )- As 
discussed in Ref. [f|, the slow running of the coupling 
suggests extracting the beta function at each (/3,K C ) from 
a linear fit of the inverse coupling 



u{L) = c + ci log 



L 



(5) 



With this parametrization, Co gives the inverse coupling 
u(L = 8a), while c\ is an estimate for the beta function 
(3 at this coupling. 

For a first look, we fit the data points for all L to 
extract the slopes at the given bare parameters, ignoring 
discretization errors that must be inherent in the smallest 
lattices. These fits are shown in Fig. Q] Values of the 
beta function /3(u) obtained from these fits are plotted 
as a function of u(L = 8a) in Fig. [5J One can see that 
the results for (3q = and for (3q = 0.5 are consistent 
with each other. Also shown are the one- and two-loop 
approximations from the expansion 



p(u) 



2bi 
16tt 2 



26 2 



1 



(16tt 



2^2 



(6) 



where h = 13/3 and b 2 = -194/3. 

The assumption behind the linear fits is that j3 is small 
so that u(L/a) changes very slowly with the volume; this 
behavior is apparent in Fig. [T] Nonetheless, where there 
are four or more volumes to be fitted for a given coupling 
the x 2 is n °t particularly good. Corrections to the ap- 
proximate Eq. ([5]) come from discretization errors, as well 



as from the slight deviation from constancy of the contin- 
uum beta function over the range of volumes. We defer 
consideration of discretization errors to Sec.fV] Variation 
of the continuum beta function gives rise to higher pow- 
ers of logL. Adding the next-to-leading term, at each 
bare coupling we fit 



u(L) — c + ci logL/8a + c 2 (logL/8a)^ 



(7) 



From the definition of a beta function it follows that c\ 
continues to provide an estimate for j3 at u — 1/ g 2 (L = 
8a). The results of these fits are shown as empty symbols 
in Fig. [21 It is evident that there is only a small change 
compared to the linear fits of Eq. ([5]). In our analysis of 
discretization error, given below, we use only the single 
logarithm. 

At the two strongest couplings the beta function plot- 
ted in Fig. [3] is positive, indicating the existence of an 
infrared fixed point. As we shall see, this conclusion does 
not survive consideration of the discretization error. It is 
possible that the beta function remains negative even at 
the strongest couplings studied. Hence, we cannot rule 
out a "walking" scenario wherein the beta function comes 
close to zero, but never actually turns positive. 

Nevertheless, Figs. [T] and [2] demonstrate an important 
qualitative result, which will be unaltered by any anal- 
ysis of lattice error: The coupling constant runs more 
slowly than one-loop perturbation theory over the entire 
range of bare parameters. This is in marked contrast to 
QCD, where the running coupling runs faster than the 
one-loop (or the two-loop) beta function as one enters 
strong coupling. 
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FIG. 2: Left: Beta function ${u) plotted as a function of u(L = 8a). The squares are from the fie — data while the circles 
are from /9e = 0.5. Results are extracted from the linear fits l[5]l. as shown in Fig. [T] Plotted curves are the one-loop (dotted 
line) and two-loop (dashed line) beta functions. Right: Filled symbols as on left; empty symbols derive from fits to Eq. l[7[l. in 
which a log 2 term has been added. The empty symbols have been slightly displaced to the right. No correction has been made 
for discretization errors. 



IV. MASS ANOMALOUS DIMENSION 

We derive the mass anomalous dimension from the 
scaling with L of the pseudoscalar renormalization factor 
Zp. The latter is calculated by taking the ratio 



0.3 



fp(L/2Y 



(8) 



fp is the propagator from a wall source at the t = 
boundary to a point pseudoscalar operator at time L/2. 
The normalization of the wall source is removed by the 
/i factor, which is a boundary-to-boundary correlator. 
The constant c, which is an arbitrary normalization, is 
1/V2 in our convention. 

We present in Tables [V] and I VII the results of calculat- 
ing Zp in our runs with /3 6 = 0.5 and with /3 6 = 0; we 
plot them in Fig. [3] 

The slow running suggests [f| that we may extract 7 m 
by applying the approximate scaling formula 



Z P (L) = Z P (L ) 



(9) 



that is, from the slopes of the lines drawn in Fig.[3J These 
linear fits are analogous to Eq. ([5]): 



log Z P (L) = c +cilog- 



(10) 



Again, we begin by fitting to all volumes at a given bare 
coupling simultaneously, leaving for later the discretiza- 
tion errors due to the smallest lattices. The individual 



0.2 




0.1 



FIG. 3: The pseudoscalar renormalization constant Zp 
vs. L/a. The crosses are from simulations with fie = 0, 
/3 = 5.8 to 4.4. The circles are from simulations with f3e = 0.5 
(TableE): top to bottom, /3 = 3.5, 2.5, and 2.0. The straight 
lines are linear fits to each set of points at given (/3, /?e); the 
slope gives — y m . The hypothetical dotted line corresponds to 



7r. 



= 1. 
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TABLE V: Pseudoscalar renormalization factor Zp evaluated at the couplings (/?, k c ), with j3g = 0.5, for lattice sizes L. 



p 








Z P 






L — 6a 


L = 8a 


L = 10a 


L = 12a 


L = 14a 


L = 16a 


3.5 


0.2171(2) 


0.1923(8) 


0.1789(6) 


0.1680(5) 






2.5 


0.1787(4) 


0.1560(5) 


0.1409(6) 


0.1325(6) 


0.1263(10) 


0.1192(8) 


2.0 


0.1579(6) 


0.1371(6) 


0.1245(10) 


0.1183(10) 




0.1024(10) 



TABLE VI: Pseudoscalar renormalization factor Zp for the strong-coupling cases with j3e ~ 0. Data for L — 10a, 14a are from 
the new simulations listed in Table [Tit while the other data are from Ref. [J]. 



p 








Z P 






L — 6a 


L = 8a 


L = 10a 


L = 12a 


L = 14a 


L = 16a 


4.8 


0.2246(23) 


0.1981(15) 


0.1824(9) 


0.1716(10) 






4.6 


0.2127(14) 


0.1808(16) 




0.1518(14) 




0.1340(6) 


4.4 


0.1888(18) 


0.1631(16) 


0.1447(6) 


0.1311(13) 


0.1256(9) 


0.1163(13) 



fits are shown as straight lines in Fig. [3J The outstanding 
feature of Fig. [3] is that in no case does the slope of any 
data set approach —1, meaning that 7 m (Fig. [4]) never 
reaches unity. This qualitative observation survives all 
further analysis. 

Allowing for running of the coupling, we have also con- 
sidered a fit function analogous to Eq. 0, 



log Z P (L) = c + ci log L/8a + c 2 (log L/8a) 



(11) 



For both fits the mass anomalous dimension at g 2 (L = 
8a) is given by — c\. We show the results of both fit types 
in Fig. 01 plotted against the running coupling g 2 (L = 
8a). It is apparent that the result for 7m, (g 2 ) changes 
little when the fit is broadened, or in other words, that 
scaling violations due to the nonzero beta function are 
small. 

A comparison of results for the two lattice actions 
shows that there is some disagreement. The two 
strongest-coupling points obtained with (3$ — lie well 
above the line connecting the (3 e = 0.5 points. As we 
shall see, analysis of discretization errors reduces this dis- 
crepancy. Moreover, the rough bound 7 m ^0.5 will be 
strengthened by extrapolation to the continuum limit. 



V. ESTIMATING DISCRETIZATION ERROR 

Before showing our extrapolations to the contin- 
uum/infinite volume limit, let us make some general com- 
ments. 

Lattice studies of quantum chromodynamics approach 
the continuum limit in the weak-coupling regime. QCD 
possesses a physical scale A. A systematic approach to 
discretization error gives an expansion of physical quan- 
tities in powers of Aa, with which one may extrapolate 
to the a —¥ limit. Perturbation theory gives a guide 



to the coefficients in the expansion. The L — > oo limit 
is taken separately; it is governed by the mass gap M of 
the theory and the corrections are functions of ML. 

The theory considered here differs from lattice QCD 
in several essential features. In the first place, the cou- 
pling runs much more slowly than in QCD. After all, the 
theory was selected for study because this is true in the 
two-loop beta function; our numerical results bear this 
out nonperturbatively. In the second place, we seek the 
putative infrared fixed point in strong coupling. Along 
with a strong SF coupling at the volume scale L, then, 
we have slow running to a still-strong coupling at the 
lattice scale a. There is a perturbative scale Ai that 
emerges from one-loop running, but it never comes into 
play. Much as in truly conformal theories, discretization 
errors are functions of a/L — they are indistinguishable 
from finite- volume corrections. It is clear from this that 
any perturbative estimate of discretization error is point- 
less. 

Nonetheless, we offer a discussion of the fermion contri- 
bution to the one-loop SF coupling in the appendix. We 
point out that, even were we to work at couplings where 
one-loop perturbation theory is valid, for our lattices the 
0(a 2 /L 2 ) corrections are so large that they cannot be 
disentangled from 0(a/L). This means that any fit func- 
tion used to extrapolate to a/L = constitutes a model 
rather than an expansion about the limit. 

Since our main result is the upper bound on the anoma- 
lous dimension 7 m , we begin with the analysis of dis- 
cretization error in this quantity. We focus on the strong- 
coupling data displayed in Tables Ivl and I VII The /3e = 0.5 
data show the dramatic departure of 7(ff 2 ) from the one- 
loop line (see Fig. Comparison to the fle = shows 
an apparent discrepancy between the two actions that 
begs explanation. We present several extrapolations to 
a/L = 0, which result in a systematic shift downwards of 
the results. At the same time, the results from the two 
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FIG. 4: Left: Mass anomalous dimension "fig 2 ) plotted against g 2 (L — 8a). The squares are from the /3a = data while the 
circles are from fie = 0.5. Results are from the linear fits shown in Fig. [3] The line is the one-loop result. Right: Filled symbols 
as on left; empty symbols derive from fits to Eq. in which a log 2 term has been added. The empty symbols have been 

slightly displaced horizontally. No correction has been made for discretization errors. 



actions are brought closer to agreement, both through 
the propagation of statistical error in the extrapolation 
and through uncertainty surrounding the method of ex- 
trapolation to a/L = 0. 

Similar analyses of the beta function are inconclusive. 
The dependence of the running coupling g 2 on lattice size 
L is irregular and not amenable to smooth extrapolation 
to the continuum. In the end, we have to let our lattice 
results stand as they are. 



A. Anomalous dimension 



The general behavior we expect for Zp(L) is 



Z P (L) = Aexp 



dt 



Ira (g 2 (tL)) 



Pi- 



(12) 



For a slowly-running coupling, the first term simplifies to 



A(L/a) 



-7 m (9 2 ) 



The original procedure [35j is to extract 



the step scaling function ap(g 2 ,s) from the ratio of Zp 
measured on two lattices, 



a P (v, s) 



Zp(sL) 



9 2 (L)=v 



Z P {L) 

Absent scaling violations, this is just 
a P {g 2 ,s) = s-^^ 



(13) 



(14) 



The traditional analysis 3 continues by computing 
ap(g 2 ,s) at fixed s for several values of L and extrap- 
olating L to infinity. Let us begin our analysis by doing 
variations on the traditional fit. 



1. Extrapolation of the step scaling function 

At 4 P = 4.4, 4.6, and 2.0, 2.5, we can do the most 
traditional fit of all: fix the ratio s = 2 and compare ap 
for the pair L — (6a, 12a) with L = (8a, 16a). The results 
for 



7m (L) 



log gp 
log 2 



(15) 



are plotted in Fig. [5] against L of the smaller lattice. 

In all four cases, there is a discernible L dependence in 
the data. Now we face the problem of how to extrapolate 
in L. As noted above, we are away from any perturbative 
limit, and so theory cannot guide us. Choosing models, 
we perform extrapolations linearly in L, 



7 m (£) = 7™ + C— , 



(16) 



3 The truly traditional analysis keeps g 2 fixed by varying the bare 
couplings as L is changed. As we do throughout this paper, we 
take advantage of the slow running to do our analyses at fixed 
(/3,/9g). The variation of g 2 is negligible. 

4 We sometimes refer to the data sets by their f3 values. The tables 
show that there is no confusion if we omit mention of j3g . 
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(a/L) 



FIG. 5: 7 m from "traditional" fits to pairs of lattices (L,2L) = (6a, f 2a) and (8a, I6a), plotted versus a/L (left) and versus 
(a/L) 2 (right). The extrapolated values are shown near the origin. Couplings (/3,/3e) are as shown. 



or quadratically, 



7m(£)=7m + C(|y 



(17) 



The fits, with their extrapolations to a/L = 0, are shown 
in the figure and tabulated in Table IVlTl Because our fits 
are of two points to two parameters, there is no goodness- 
of-fit criterion to invoke: The data cannot express a pref- 
erence between the models. 

The extrapolations to a/L = show a (mostly) down- 
ward trend compared to any individual ap. The result- 
ing error bars are of course larger but 7 m never exceeds 
0.5 for any of the couplings examined. We defer further 
discussion of this sort until we reach the results of our 
all- volume analysis below. 

In principle, we can judge the quality of the fits if we 
add data points, that is, more lattice volumes. To this 
end we have run simulations for lattices of size L = 10a 
and 14a at two of the couplings, j3 — 2.5 and 4.4. This 
gives us three independent pairs of volumes from which 
to calculate a step scaling ratio. We use the pairs L/a — 
(6, 12), (8,14), and (10,16) to calculate a P (g 2 ,s) for s = 
2, 7/4, and 8/5, respectively, and estimate 



7m (L) = 



log (Jp{g 2 ,s) 
logs 



(18) 



as a function of the smaller L of each pair. Fig. |6] shows 
the result, together with the extrapolations to infinite 
L — now two-parameter fits to three points, which are also 
tabulated in Table I VIII The quality of the fits is good in 
all cases, so (unfortunately) the data do not express a 
choice between 1/L and 1/L 2 corrections. 



p 


L pairs 


linear fit 


quadratic fit 


4.4 


(6,12),(8,16) 


0.37(10) 


0.44(5) 


4.6 


(6,12),(8,16) 


0.27(7) 


0.36(4) 


2.0 


(6,12),(8,16) 


0.43(7) 


0.43(4) 


2.5 


(6,12),(8,16) 


0.26(5) 


0.33(3) 


4.4 


(6,12), (8,14), (10,16) 


0.35(6) 


0.42(4) 


2.5 


(6,12), (8,14), (10,16) 


0.23(4) 


0.31(3) 



TABLE VII: "Traditional" extrapolations of j„ 
pairwise step scaling function ap. 



from the 



2. All-volume extrapolations 

Now we present a procedure for using data from all 
volumes simultaneously to arrive at the extrapolated 7 m 
at each value of the bare coupling. The most straight- 
forward approach would be to take the measurements of 
Zp and fit them to 

log Z P {L) = A + 7m log(a/L) + C j (19) 

or to 

log Zp{L) =A + lm log(<z/L) + c (|) 2 . (20) 

We find, however, that our data do not distinguish well 
among the basis functions in these fits and hence the de- 
terminations of the parameters are unstable. We there- 
fore adopt an alternative approach to estimation of the 
discretization error. 
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FIG. 6: 7 m from fits to three pairs of lattices (L, L') = (6a, 12a), (8a, 14a), and (10a, 16a), plotted versus a/L (left) and versus 
(a/L) 2 (right). The extrapolated values are shown at the origin. 



For each bare coupling ft we have measured Zp on a 
set of volumes L\ < Li < ... < Lpf. The comparison 
in Fig. [4] shows that the running of the coupling is slow 
enough to be neglected over the range L\<L< L^. The 
only obstruction to extracting 7 from a fit to Eq. (j!0[) is 
then discretization error. Since this error is most pro- 
nounced at the smallest L's, we can improve our esti- 
mate of the true 7 by successively dropping the smaller 
volumes. This observation is the basis for the following 
continuum extrapolation. 

Let Cq™ , be the parameters of the linear fit (|T0|) 
where only the largest volumes L ni L n+ i, . . . , Ln are re- 
tained. The slope is an estimate for —7 which 
we take to be a function of a/L n . For example, if 
{L\, . . . , £4} = {6a, 8a, 10a, 12a}, a fit using all vol- 
umes gives an estimate of 7 for a/L — 1/6; drop- 

ping the smallest volume gives the next estimate c\ for 
a/L = 1/8; and dropping the two smallest volumes gives 

the estimate c\ for a/L = 1/10. Finally we extrapo- 
late these estimates to the continuum limit, a/L — > 0, by 
fitting to either Eq. ([16]) or Eq. (fl7)) . 

When performing the continuum extrapolation, we 
must take into account that the results of the suc- 
cessive fits are correlated. For the example above, with a 
maximum of two volumes dropped, it can be shown that 

the covariance matrix C mn — cov ^Cj^Cj" 1 ^ is 

/AW AW AW \ 
C = AW AW A( 2 ) , (21) 
\ AW AW A<3) / 



where A^™) = a 2 f c i J i s the variance of each q™ . 

We show in Fig. [7] the linear and quadratic continuum 
extrapolations for ft = 2.5 and 4.4, where we can work 
from six volumes. In each case the fit has very low x 2 
as long as we exclude the final (leftmost) data point, 
which is based on L/a = 14 and 16 only. Since its error 
bar is large, omitting this point from the fit leaves the 
extrapolation unchanged. This is a feature of all the fits 
in the all- volume extrapolations. 

To give an idea of the stability of our fit procedures, we 
compare them all (as available) in Fig. [51 We display the 
results of the all-volume extrapolations for all bare cou- 
plings in Table I Villi and we plot them in Fig. [HI Every 
point lies lower than the corresponding point in Fig. 21 
which shows the result of the naive analysis. Correc- 
tion for discretization error has pushed our bound on j m 
downward. In addition, the gap between the two actions 
(fte = vs. 0.5) has narrowed, for two reasons. One is 
the larger error bars compared to Fig. [4j the other is the 
systematic uncertainty hovering over the choice between 
linear and quadratic continuum extrapolations. 



B. Beta function 

We would like to perform a similar analysis for the 
running coupling. Unfortunately, the data do not per- 
mit the extraction of extrapolated results carrying useful 
uncertainties. The reason why can quickly be seen from 
Fig. [TDJ The data at each individual bare coupling are 
irregular functions of L. This may mean that the error 
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FIG. 7: Examples of linear and quadratic continuum extrapolations of the mass anomalous dimension "fig 2 ) derived from 
all-volume extrapolations. Values of (j3,/3e) are as shown. 



p 


L values 


linear fit 


7m 

quadratic fit 


4.8 


6,8,10,12 


0.27(7) 


0.32(4) 


4.6 


6,8,12,16 


0.35(4) 


0.40(3) 


4.4 


6,8,10,12,14,16 


0.39(4) 


0.43(3) 


3.5 


6,8,10,12 


0.24(5) 


0.29(3) 


2.5 


6,8,10,12,14,16 


0.26(3) 


0.33(2) 


2.0 


6,8,10,12,16 


0.33(4) 


0.37(3) 



TABLE VIII: Continuum limit of 7™ from all- volume extrap- 
olations. 



bars for the individual data points are underestimated. 
The SF observable generally has a very long autocor- 
relation time. For lattice sizes L > 10a, we determined 
an error bar by fitting four runs to a common constant, 
and collecting data until the x 2 P er degree of freedom 
from this fit fell below about two. This may not have 
been good enough. From a practical point of view, the 
pattern of our data makes analysis of the discretization 
errors quite difficult. 

We shall treat only two data points: /3 = 2.0 and 2.5, 
with /?6 = 0.5. These are our strongest couplings, and as 
seen in Fig. [2] they hint at a positive beta function and 
hence an infrared fixed point. 

At /3 = 2.0 we have four volumes: L = 6a, 8a, 10a, 12a. 
(We have discarded L = 16a as discussed above.) For 
an almost traditional analysis we take them as two pairs, 
(6a, 10a) and (8a, 12a), and for each pair we calculate the 



discrete beta function (DBF), 
1 1 



B(u,s) 



9 2 {sL) g 2 {L)' 



9 2 {L) 



(22) 



The scale factor s is 5/3 for the first pair and 3/2 for 
the second pair. The two DBFs can be combined if we 
rescale them according to 



R(u, s) 



B{u,s) 
logs 



(23) 



This rescaled DBF approximates the usual beta function 
(3(1/ g 2 ) when the running is slow @. We can plot it 
against the smaller L in each pair and extrapolate to 
a/L = either linearly or quadratically, much as we did 
for 7 m above. 

At f3 = 2.5 we have six volumes. Following our pro- 
cedure for 7 m , we can first calculate the DBF for s = 2, 
using the pairs (L,2L) — (6a, 12a) and (8a, 16a). Ex- 
trapolations to a/L = then follow. Additionally, we 
calculate the rescaled DBF R(u, s) for the three pairs 
(L,L') = (6a, 12a), (8a, 14a), and (10a, 16a), and ex- 
trapolate it to a/L=0. Unfortunately, the irregularity 
of the data seen in Fig. [TU] renders the three- pair fit use- 
less: Trying to fit a straight line through the three points 
gives an enormous x 2 ■ The two-pair fits, of course, have 
no degrees of freedom and hence no x 2 ■ 

At both couplings we can do all- volume extrapolations 
as for Zp, mutatis mutandis. We fit to Eq. ([5]) to ob- 
tain coefficients , c^™' , where n indicates the smallest 
volume L n retained in the fit. Finally we extrapolate to 
a/L = either linearly or quadratically. 
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FIG. 9: Continuum limit of the mass anomalous dimension 
y{g 2 ) derived from all- volume extrapolations (see Table |VIH1 
and compare Fig. [4|. The squares are from the /3e — data 
while the circles are from f3e = 0.5. Empty symbols assume 
discretization errors proportional to a/L while filled symbols 
assume (a/L) 2 . The filled symbols have been slightly dis- 



placed horizontally. 



FIG. 8: Comparison of the results of all continuum extrapola- 
tions for the anomalous dimension at four bare couplings. For 
each coupling, we plot in sequence, top to bottom: two-pair 
extrapolations — linear and quadratic (Fig. [5}; three-pair ex- 
trapolations (where available) — linear and quadratic (Fig. [6]); 
all-volume extrapolations — linear and quadratic (Table IVIIII 
and Fig. . 



For j3 = 2.0 the all- volume extrapolations give rea- 
sonable x 2 an d improve on the pairwise extrapolations, 
as may be seen in Fig. QTJ Nonetheless, the error bar for 
each extrapolation has grown by a factor of 4-6 compared 
to the raw fit shown in Figs. [T] and [21 which wipes out 
any sign of a zero crossing. For f3 — 2.5 the all-volume 
extrapolations give positive results but, again, the error 
bar has grown by a factor of 4-6 and the central values 
lie barely 1<7 above zero. We therefore cannot confirm 
the existence of an infrared fixed point. It is quite pos- 
sible that the beta function approaches zero as seen in 
Fig. [5] but then veers away, staying negative and lead- 
ing to strong coupling at large distance. This would be 
exactly the behavior conjectured for walking technicolor. 




FIG. 10: A blown-up view of Fig. [T] showing the strong- 
coupling portion of our data sets for the running coupling 
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FIG. 11: Results of continuum extrapolations for the beta 
function at the two strongest couplings. For each (/3,f3e) we 
plot in sequence, top to bottom, the two-pair extrapolations — 
linear and quadratic — and the all-volume extrapolations — 
linear and quadratic. The vertical dotted line is the one-loop 
value. 



VI. DISCUSSION 

In this paper we have continued our study of the SU(3) 
gauge theory with two flavors of sextet fermions. A new 
term in the lattice action allowed us to explore a much 
wider range of the renormalized coupling. The old and 
new actions give rise to consistent results for the beta 
function where they overlap. For the mass anomalous 
dimension there are some disagreements; the source of 
the disagreement must be the proximity of the first-order 
phase transitions — a lattice artifact — in the old action. 
The discrepancies are reduced when we extrapolate to 
the continuum limit. 



A. An infrared fixed point? 

The question of the existence of an IRFP in this model 
has been addressed in the literature in several ways. A 
first hint comes from qualitative features of the phase 
diagram of the lattice theory [3, 0] ■ I n QCD on a finite 
lattice, there is a confining phase for sufficiently strong 
coupling within which one may tune k — > k c in order 
to obtain a massless pion. Then one can increase the 
size of the lattice while tuning towards weak coupling on 
the K c {0) curve, thus approaching the continuum limit 



within the confining phase. In the sextet theory, on the 
other hand, the confining phase is bounded by a first- 
order phase transition that does not permit definition of 
k c . It is not at all clear that a continuum limit in the 
confined phase would yield a theory with finite masses, 
let alone a chiral, massless limit. The only alternative 
might be to begin in the weak-coupling, non-confining 
phase, where one could take a non-confining continuum 
limit. 

This feature of the phase diagram was first seen in 
QCD with a large number of (color-triplet) flavors by 
Iwasaki et al. [42|, El], and explored more recently in 
Ref. [3]. This may be a sign of entry into the con- 
formal window in these theories. Turning to staggered 
fermions, the authors of Refs. [45|,[46[ have recently noted 
the absence of spontaneous breaking of chiral symmetry 
in strong-coupling (triplet) QCD when the number of fla- 
vors is large. 

Even if one ignores the issue of a light pion, one can ask 
whether any kind of continuum limit can be taken in the 
confining phase. On a finite lattice, the first-order phase 
boundary at strong coupling connects to the boundary 
between strong-coupling and weak-coupling phases that 
may be interpreted as the finite-temperature phase tran- 
sition. In a confining theory, this transition must move 
towards weak coupling as the lattice is enlarged, in such 
a way as to give a finite transition temperature in the 
continuum limit. In the theory with thin-link Wilson 
fermions, we found }2| that this motion stalls near the 
weak-coupling n c curve, as if the finite-temperature tran- 
sition is avoiding the basin of attraction of an IRFP. More 
detailed study of this question in the fat-link theory 0] 
revealed a slow motion of the transition with lattice size; 
one would need to study larger lattices to see whether this 
motion eventually gives correct scaling in the continuum 
limit. In fact, because the perturbative running of the 
coupling is so slow [47j, it will take enormous lattices to 
settle the matter. Kogut and Sinclair [48U5l| have been 
studying the same question with staggered fermions; the 
results are, so far, similarly inconclusive. 

Fodor et al. [521 - 156} have studied the staggered- fermion 
theory with extensive simulations on large volumes. On 
each lattice, they use quark masses that are large enough 
that the volumes are effectively infinite. They then test 
alternative scaling hypotheses for the mass spectrum, the 
chiral condensate, and the string tension as a function 
of the quark mass. Their analysis favors confinement over 
conformal physics, meaning a beta function without an 
IRFP, which might behave as a walking theory. 

As we have seen, our Schrodinger-functional results are 
not precise enough to nail down a fixed point. Statistical 
fluctuations preclude unambiguous extrapolation to the 
continuum limit. All we can say is that the beta function 
is smaller than the one-loop value, so that the theory runs 
slowly. 
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B. The mass anomalous dimension 

As we have seen, j m first follows the one-loop curve 
and then at stronger couplings it levels off. This is in 
line with what we found in the SU(2) and SU(4) theo- 
ries @,@|- Our simple fits, without taking lattice error 
into account, give a bound 

7m £ 0.45 . (24) 

At weak couplings we did not carry out a continuum 
extrapolation. At our strongest couplings, our analysis 
of discretization error leads to larger error bars than our 
first analysis, but also to a general movement downwards 
so that our final result is the same bound ([24)) . 

The issue of scheme-dependence always arises at this 
point. The value of j m at a fixed point, if there is one, 
is scheme-independent. We can say more than this, how- 
ever. A bound like Eq. (j2"4"]l is evidently invariant under 
any redefinition g — > g'{g), that is, it is entirely scheme- 
independent. 

A successful model of walking (extended) technicolor 
must have a slowly-varying coupling, but no infrared 
fixed point; it must also have a large mass anomalous 
dimension, j m ~ 1. While our results do not rule out 
walking in the sextet theory, the smallness of 7 m makes 
this theory unsuitable for walking technicolor. 
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Appendix: One loop analysis 

In the one-loop approximation the SF coupling, defined 
in Eq. ([3]), takes the form 



Here go is the bare coupling from Eq. ([2]), and the one- 
loop contributions are Eg from the gauge and ghost 
fields 5 and Ef from each flavor of the sextet fermions. 
Following the definitions 6 and methodology of Ref. [2(| , 
two groups [HrHH have recently calculated the fermion 
piece Ef for a number of gauge groups and representa- 
tions and for a variety of boundary conditions and back- 
ground fields (as imposed in the SF method). On general 
grounds, one can expand 

Ejr(£) = r + s log(L/a)+ri— + si—log(L/a) 

2 2 

+r 2 ^2 + s 2 lo &( L / a )j2 +■■■ ■ ( A - 2 ) 

The continuum limit consists of the first two terms on 
the right-hand side. The coefficient sq of the logarithmic 
term is the universal contribution of a fermion field to the 
one-loop beta function. The constant r$ depends on the 
choice of the background field, i.e., it introduces scheme 
dependence. 

All other terms on the right-hand side of Eq. (|A.2[) 
constitute the discretization error. In the presence of a 
clover term with csw = 1 we have s% = 0. The remaining 
terms do not vanish in general. The linear discretization 
error (rj term) can be cancelled by adding a boundary 
counterterm for the gauge field. As for the quadratic 
discretization error (the r 2 term), one can suppress it by 
a judicious choice of the background field |59l - [61l |. The 
outcome of adopting both measures is an improved one- 
loop behavior, where s% = r\ = r 2 = 0. 

The above strategy is natural for QCD, where the bare 
coupling is such that physics at the cutoff scale is typ- 
ically well inside the perturbative regime. Here we are 
faced with a drastically different situation. Because of 
the slow running, in order to reach the interesting strong 
values of the renormalized coupling, the bare coupling 
must already be strong. As a result, we can no longer rely 
on perturbation theory at the quantitative level. Suffice 
it to mention that, in the vicinity of a two-loop fixed 
point, the two- loop contributions must be as important 
as the one-loop! Besides, the main thrust of the present 
paper is the calculation of the mass anomalous dimen- 
sion for sextet fermions, for which no perturbative study 
of discretization errors has been carried out. 

The conclusion, obviously, is not that discretization 
errors are to be ignored; as in any lattice simulation, 
a continuum extrapolation is mandatory. The point is 
that the continuum extrapolation must be done non- 
perturbatively. This is what we have attempted, with 
mixed success, in Sec. |V] 

We now proceed to a discussion of Ef for sextet 
fermions, in order to get some idea of the magnitude of 
this piece of the discretization error in the weak-coupling 



(A.l) 



T,q(L) has not been calculated for nHYP links. 
6 Our E F (L) is called p M in Ref. [H. 
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regime. We have computed Ef by following closely Ap- 
pendix A of Ref. [2(J. Inputs to the calculation are the 
fermion twist phase, set to 8 = ir/b, and the classical 
background field, which, in turn, is determined by the 
SF boundary data. (nHYP smearing leaves the classi- 
cal background field unchanged.) The spatial links on 
each time boundary are specified by three abelian phases 
4>i = 'Piiv): i = 1,2,3, which satisfy the constraint 
4>i + 02 + 03 — 0. The corresponding six phases of the 
sextet representation are given by <pij = <pi + <fij, with 
1 < i < j < 3. Ef is written as a sum over the lattice 3- 
momentum p. For each p, the fermions' contribution can 

be represented as d/drjtr logM , where M is a two- 

77=0 

by-two matrix that depends on p and on the boundary 
conditions (2p| ]. 

We plot Ef as a function of L in Fig. [T3J For 
comparison, we also show the continuum result Ef = 
r +b/(12n 2 ) logL, where the scheme-dependent constant 
ro has been extracted from a fit of Ef to the six terms 
shown in Eq. (|A.2[) . In this fit, the coefficient so is fixed 
at the continuum value and s± is fixed to zero (see above). 
The other coefficients are listed in Table IIXI 

In order to get an idea of the magnitude of discretiza- 
tion errors in the beta function, we plot Ef against 
log(a/L) in the second panel of Fig. |T2} here we limit 
the range of L to 6a < L < 16a to match our numer- 
ical simulations. This plot is analogous to Fig. [1] (but 
note that here we display only the contribution of a sin- 
gle sextet fermion). The continuum result now shows 
up as a straight line, whereas the lattice result exhibits 
deviations from linearity. It can be seen that much of 
the discrepancy between the continuum and lattice re- 



sults can be accounted for by an additive constant. As 
noted above, a constant shift ro does not represent a dis- 
cretization error, but, rather, scheme dependence. The 
discretization error in the beta function will appear as 
a difference in the slopes of the two lines. Using the 
method of least squares to estimate the average slope of 
the lattice data we find the value 0.0446, which amounts 
to a 5% deviation from the correct continuum value 
5/(12tt 2 ) ~ 0.0422. Multiplying by two for the number 
of flavors, we conclude that the fermions' contribution to 
the discretization error in the beta function is in abso- 
lute terms roughly 0.005. A glance at Fig. [5] shows that 
this is a less than half the statistical error plotted for the 
weak-coupling data points. 

The smallness of the deviation in the slope is connected 
to the coefficients r\ and r-i of the linear and quadratic 



TABLE IX: Free coefficients in the fit JA2 



ro 


0.1569 


ri 


-0.1914 


ri 


0.4169 


S2 


0.0086 



lattice corrections. The two coefficients turn out to have 
opposite sign and magnitudes that cause their contribu- 
tions to the slope to cancel. Thus one cannot really jus- 
tify an expansion in powers of a/L in extrapolating to 
the continuum limit from our data. As we have stressed 
above, any given extrapolation is a model. 
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